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SUMMARY 

Tuning maps are an aid in the controller tuning process because they 
provide a convenient way for the plant operator to determine the consequences 
of adjusting different controller parameters* In this application the maps 
provide a graphical representation of the effect of varying the gains in the 
state feedback matrix on startup and load disturbance transients for a three 
capacity process. Nominally, the three tank system, represented in diagonal 
form, has a Proportional-Integral control on each loop. Cross-coupling is 
then introduced between the loops by using nonzero off-diagonal proportional 
parameters. Changes in transient behavior due to setpoint and load changes 
are examined by varying the gains of the cross-coupling terms. 


INTRODUCTION 

Tuning of control systems to achieve desired transient response charac- 
teristics to setpoint and load changes has long been an important issue. In 
1942, Ziegler and Nichols proposed the first structured approach to the tuning 
of Proportional-Integral-Derivative (PID) controllers to appear in the litera- 
ture (ref. 1). This technique presented a simple way to achieve one-quarter 
damping in response to load disturbance upsets in industrial-type single- 
input-single-output processes. Since then, many forms of the algorithm have 
appeared with the emphasis on different types of responses (refs. 2 to 6). As 
modern control techniques have become more popular in practice, some tuning 
methods for multivariable systems have begun to appear in the literature 
(refs. 7 to 9 ) . 

The original, pre-Ziegler-Nichols tuning method, sometimes known as 
tweaking or trimming, is still very much in use in industry. It consists of 
an operator observing a transient, adjusting the controller gains based on 
experience, and rerunning the transient. This method is often limited in 
practice by the infeasibility of running exploratory step responses on-line 
and by the length of the dominant time constant of the plant which may be 
hours for some processes. During normal plant operation, an operator is 
generally more concerned with load disturbance rejection than with response to 
setpoint changes. However, since it is often impractical to generate load 


Summer High School Apprenticeship Research Program at NASA Lewis 
Research Center. 



upsets for testing purposes, tuning of a controller to compensate for loads 
may be difficult. Thus, operator dissatisfaction with closed loop performance 
is common and results in about half of all industrial control loops being run 
on manual (in open-loop mode) (ref. 10). 

Operator aids known as tuning maps were presented in 1962 (refs. 11 
and 12). These maps consisted of transient responses of a specified shape 
(one-quarter damping, for instance) plotted in the plane comprised of the nor- 
malized integral and derivative gains. The loop was closed around a multiple 
lag plus deadtime process with a PID controller. The integral and derivative 
terms were fixed and the proportional value was varied until the step response 
had the desired characteristics. The transient was displayed at the appropri- 
ate point in the normalized Integral-Derivative plane and the proportional 
value which produced the response was shown alongside it. The map allowed the 
operator of a process with the designated transfer function to choose the PID 
values directly from the map of the desired response curve shape. 

The advantage of the tuning map over the Ziegler-Nichols-type tuning 
algorithms is that it gives a large choice of PID par am eters rather than a 
single set for the desired response shape. Since the Ziegler-Nichols algo- 
rithm seeks to achieve only one feature of the response (damping) as its 
objective, other parameters such as period of oscillation, percent overshoot, 
etc., are not considered. Operators who must meet several of these con- 
straints simultaneously need the flexibility to choose initial controller 
values which will come close to attaining their goals at once. 

An extension of the tuning map idea which offers even more flexibility is 
the three-dimensional tuning map. This map again contains transient plots at 
specific points in the normalized Integral-Derivative plane but with many pro- 
portional values used for each Integral-Derivative pair. Thus, the map gives 
the additional information of the direction of adjustment for the tuning 
parameters to attain the desired response. This technique allows an overall 
response curve to be chosen rather than simply a single feature (ref, 13). 

All previous work has been applied to single-loop processes using class- 
ical control methodologies. Our work involves using tuning maps for a state 
feedback controller to demonstrate the effect of introducing interaction 
between multiple loops of a formerly decoupled process. Generally it is 
desirable to have noninteracting loops for ease of tuning, but, if the inter- 
action provides the potential for improved system response and the method for 
tuning the controller is provided, the addition of interaction can be benefi- 
cial. The sample process used consisted of three noninteracting cascaded 
capacities. This system was selected arbitrarily as a representative indus- 
trial process. The scope of the work was limited to having at most one cross- 
coupling term in each multivariable controller examined. It was felt that 
exploring more complicated interactions would obscure the purpose of the maps 
by creating an unwieldy number of variables which would overwhelm rather than 
aid the operator (ref, 14). 


The remainder of the paper consists of three parts: a description of the 

process and the artificial interactions imposed on it; a summary and analysis 
of the tuning maps obtained for the process; and conclusions about the utility 
of the maps. 
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THREE CAPACITY PROCESS 


The three tank system shown in figure 1 is an example of the noninteract- 
ing three capacity process of interest. We will assume that each tank is tall 
enough with enough liquid that transient responses will not cause any tank to 
empty or overflow, and that the valve determining the flow out of the tank 
will never fully open or fully close. These restrictions will assure linear 
behavior. We will also assume that the level in each tank is regulated about 
its own setpoint. Finally, we are concerned only with the level of the bottom 
tank. The levels of the upper two are measured as intermediate variables for 
control purposes only, they are not system outputs. 

The differential equation representing the liquid level, h, inside each 
tank is a function of the inflow, q , and outflow, q out * namely, 

Ah(t) = q in (t) - q out (t) 

where A is the cross-sectional area of the tank which is assumed to be 
unity . 

We assume that the outflow is regulated by a valve which can convert a 
flow command to an appropriate opening in negligible time. Thus the transfer 
function of the valve from flow request to output is 

G v (s) = 1.0 (!) 


A Proportional-Integral controller of the form 

k r k p s + k z 

G c (s) = k p + _ = 

c s s 


( 2 ) 


is used to govern the level in the tank. The closed loop transfer function 
for any one tank is 


h(s) 


r 

_ + 


. Ar 


S 2 . . 

s + k p s + k T 


(3) 


where r represents the new setpoint for the liquid level in the tank which 
we assume is a constant, and Ar is the amount of the step change in setpoint 
(see the appendix for the derivation). Figure 2 depicts a block diagram of 
one of the tanks. The transfer function of the error in the level (actual 
minus desired) of any one tank is simply 

e(s) ■ Ar 

s 2 + k p s + k z 

assuming the disturbance input, q in (t), is unchanged. In the case of a 
setpoint change at the lowest tank only, the transfer function for the level 
in that tank is given by equation (3). As long as the capacities are 
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noninteracting, the Proportional-Integral controller for the lowest tank 
determines the transient response for the level as it is adjusted due to the 
setpoint change. 

For load disturbance upsets, Aq in , the closed loop transfer function for 
any individual tank is 


h(s) 


r 

— + 
s 


Aq in (s) 


+ k p s 


+ k T 


(4) 


where r represents the setpoint for the liquid level in the tank and we 
assume that it is a constant. The transfer function of the error in the level 
{actual minus desired) of any one tank is simply 

«< S > = — Hn< S > (5) 

s + k p s + kj 


From figure 2 it can be seen, using equations (1), (2), and (5), that 


k p s + k r 

<iout< 8 > - — A< W S > 

s 2 + k p s + kj 


(6) 


Therefore, in the case of a flow change, Av, into the highest tank, the trans- 
fer function for the level in the lowest tank is 


h(s) - _ + — 
s _2 


k p s + k : 


k p s + k 

P - — Av(s) 


s'* + k p s + k-j. b 2 + kpS + k r s 2 + k p s + k z 


which is a combination of equations (4) and (6) since the outflow of one tank 
is the flow into the next, causing the error to propagate through the three 
noninteracting tanks. 

Using the modern control technique of state feedback to, in effect, cou- 
ple the capacities, it might be possible to improve the transient response to 
setpoint or load changes by restricting or increasing the flow from the higher 
tanks to the bottom one. 
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A possible state space realization for the three tank process is 


= Ah + B q <3out + FV = A < e + r > + B q <3out + 


Fv 


Ae + Ar + B q q out + Fv = 


Ae + B q q 0 ut + Fv = 


Ch » [0 0 1J 
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where v is the unrestricted flow into the top tank. 

In order to include the integral terms of the controllers, the state vec- 
tor can be augmented by the integrator state variables, x , x , and x . 

These additional variables provide the integral of the error in the level of 
each tank. The sixth order equation is 


0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
1 0 0 0 0 0 
0 1 0 0 0 0 
0 0 1 0 0 0 


h 3 * [0 0 1 0 0 0] 


-1 

1 

0 

0 

0 
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1 

0 

0 
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0 

0 

-1 
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Sout, 


(7) 
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If we Introduce the state feedback matrix, K = [K p K j, which incorporates the 
proportional and integral gains, we can obtain a closed loop system of the 
form 
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by setting the vector q out from equation (7) to 
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k I 
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( 8 ) 


This gives us a realization of the three noninteracting capacities pictured in 
figure 1. The fact that the portions of the K matrix corresponding to the 
proportional and integral parts of the controller are diagonal keeps the tanks 
decoupled. 


Introducing interaction between the tanks in order to improve the startup 
or load disturbance rejection characteristics involves simply adding nonzero 
off-diagonal terms to the K matrix. In this work we are concerned only with 
the effect of proportional coupling on the system. The state space formula- 
tion of the interacting system is 
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h 3 = [0 0 1 0 0 0] 


(9) 


The tuning technique we will present may be understood symbolically from fig- 
ure 3. Each knob on the front of the box in figure 3 represents a propor- 
tional element of the left half of the K matrix. The tube on top provides 
the disturbance input to the uppermost tank. The setpoints for the tanks 
enter on the left. The process starts out in steady state. The output due to 
a step change in the level command or the disturbance input is a transient 
response which is plotted on a strip chart recorder to the right. Many tran- 
sients are obtained, each with the knobs adjusted differently. The plots are 
collected and presented in a way which makes it easy to see the effect of each 
tuning change on the output. 

Figures 4 to 6 show the tuning maps generated for step setpoint changes 
from the system in equation (9). Figure 7 depicts the tuning map created for 
the system in (9) due to a step load disturbance. The three-by-three grids of 
plots are arranged so that each set of graphs in each figure corresponds to 
the adjustment of the corresponding tuning knob from figure 3. In each of 
these figures, the nominal or noninteracting system (the combination of 
equations (7) and (8) with k lAi = l and k Pii =2) is used. One by one, each 
proportional element of the K matrix is varied through a range of values and 
then reset to its nominal value. Each set of graphs depicts step responses to 
setpoint or load changes. Within each graph, the corresponding proportional 
gain is increased from upper left to lower right. In figure 4 the setpoint of 
the bottom tank only is changed; in figure 5 the desired levels of the bottom 
two tanks are changed simultaneously; and in figure 6 all three requests are 
altered concurrently. 
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RESULTS 


The tuning maps in figures 4 to 6 each represent different situations. 
The case presented in figure 4 concerns a setpoint change in the bottom tank 
only. Consequently, since the process is originally in steady state, feeding 
the error from either of the upper two tanks to any other loop has no effect 
on the transient response. This is because no error has been introduced in 
the levels of the upper two tanks. The result is that the tuning maps vary 
with the element of K in the right column only. Likewise, in figure 5 where 
the command is changed in the lower two tanks, feeding the error in the level 
of the uppermost tank to either of the lower two has no effect since the error 
is zero. This results in the left column of the map being unaffected by vari- 
ations in K. Figure 6 shows the case where all setpoints are changed making 
any change in K have an effect on the step response transient. 

The tuning map in figure 7 represents the case where a step change is 
introduced in the flow into the top tank. This disturbance upBets the equi- 
librium in the tanks, and the valves controlling outflow must be adjusted to 
return the tanks to their desired levels. The plots depict the impact of 
altering the interaction between the tanks on the level in the bottom tank. 
Since the error is injected at the beginning of the flow path and its effect 
is propagated through to the end, any change in the interaction between upsets 
causes the ensuing transient responses to differ. 

One caveat concerns the settings of parameters which apparently have no 
effect on the transients in figures 4 and 5. The fact that no changes occur 

in the step response due to varying the K value does not mean that the gain 

selection is arbitrary; only when the error is exactly zero will this situa- 
tion arise. Uninformed choice of the gains could result in an unstable system 
which will manifest itself only due to noise or disturbance upsets. In this 
case, any nonzero error will cause a runaway in the unstable loop, creating a 
serious upset in the output. Since the output of this three capacity process 
is ultimately limited by flow into the top tank, the desired output steady 
state level might be achieved but the transient response will be unpredict- 
able, defeating the purpose of the tuning map. 

Also, the responses from the figures due to varying particular K ele- 

ments away from the nominal form are in general not additive. That is, alter- 
ing two or more K values at once in order to achieve, in some sense, a 
combination of two responses from the maps will not work and may, in fact, 
give an unstable result. It is entirely possible that simultaneously altering 
several gains from their nominal settings will give greatly improved results. 
However, the tuning maps presented here will not provide the information on 
how to accomplish this. 


CONCLUSIONS 

Tuning maps for a three capacity process under multivariable control 
were developed as a guide for choosing gains of a state feedback controller. 
The maps provide a concise graphical representation of the effect of altering 
a single proportional gain on startup and load disturbance transients. The 
use of the maps is straightforward and easily understood. 
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The difficulties discussed at the end of the results section are an 
indication of the complexity encountered when introducing interaction to even 
a simple system. The maps are intended to be used as a guide during the tun- 
ing process, not as a design aid. As such, they demonstrate the effect of 
simple parameter changes on step response transients and make no attempt to 
address the result of changing several gains at once. 

The determination of what constitutes the best transient response to a 
setpoint or load change basically depends upon the operator's preferences. 
There may be some constraints tied to the process's output product quality, 
but otherwise the operator is generally free to choose the response shape. 

The use of multivariable tuning maps gives the operator an extra degree of 
freedom over standard single loop tuning methods and provides the potential 
for greatly improved control. 
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APPENDIX 


The derivation of the transfer function for the transient response of 
the level in a single tank to setpoint and load changes follows. Refer to 
figure 2 for a block diagram of a single tank. 

We may begin with the differential equation for the level in the tank 

Ah(t) = q in (t) - q out (t) (Al) 


and for simplicity we assume the cross sectional area. A, of the tank is 
unity. 

If we assume that the system is in steady state before the setpoint and 
load changes, we have the relation that 

h ( t ) = 0 

which, combined with equation (Al), gives us 

= Sout^) or Sin^o") " Sout^O-) < A2 > 

where t Q - is the time just before the setpoint and load changes. 

Referring to figure 2 we see that 

<W( S ) = — — g “— (h(s) -r(s)) 

which has the particular solution 
SC 3out< S > -Sout^O-) 

W s > 


(k p s + kj. ) (h(s) - r (s)) 
Sout^o") Jss + k r 


(h<g) - r(s)) 


(A3) 


Immediately after the setpoint and load changes occur, the level in the tank, 
h(t Q +), is the same as it was just before the setpoint change. Thus it is 
clear from figure 2 that 

h(t 0 -) = h (t 0 +) = r + Ar (A4) 

where r is the new setpoint and Ar is the value of the step change in 
setpoint, equal to the error at time t Q +. 

Finally, from time t + on, we have the fact that 

r(s) = — and q in (s) = — - in ^ to - + Aq in (s) ( AS ) 

s s 

where Aq ^ is the change in load. 
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Transforming equation (Al) into the Laplace domain and substituting in 
equation (A3) gives 

sh(s) - h(t 0 -) = q ln (s) - q out (s) 

sh(s) = h(t 0 -) + q in ( s ) - - k P s + k i (h(s) _ r (s)) 

= h(t 0 -) + q in (s) - 


s s 

Soutfto") kpS + k x 


s 


k p s + kj 

h (s) + — J r(s) 


Substituting in equations (A4) and (A5) and rearranging gives 


sh(s> ♦ -ir.i. * A,, Is) - * **“ ’ r 


Qout^o") i k P S + k I 


Using equation (A2) and combining terms gives 


I h(s) = Ar ♦ q ^ (t °~ l + Aq ln (s) - SsUL £s± * 


h(s) = 


s 2 + k D s + k T 


-Ar + — - Aq in (s) + i 

s 2 + k p s + kj s 
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Figure 4.— Tuning map of the response of a threensapadly process under multivariable control to a set point change k\ the bottom tank. 
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